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Abstract 

We analyse Belle data on the decay r~ — > n~7rv r using a dispersive representation of the vector form factor which is 
consistent with chiral symmetry and preserves analyticity and unitarity exactly. We fit the unknown theoretical param- 
eters from the data, determining the values of the related low-energy observables (r 2 ^ and Cy. The implementation 
of isospin breaking effects is also discussed (TJ- 
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1. Introduction 

The vector form factor of the pion, F v (s), encodes all 
unknown strong interaction dynamics in t~ — > n~7t°v T 
decays 0. It is defined as 

< 7r-{p n -)n\ Pn0 )\dru\Q >= V2(/v -/vfW> W 

with s = (p„- + p n o) 2 . This form factor is not only rele- 
vant for the understanding of the hadronization of QCD 
currents at low energies (see [3]), but also represents a 
crucial ingredient for the evaluation of the leading order 
(LO) hadronic contribution to the anomalous magnetic 
moment of the muon, which provides a very stringent 
probe of new physics |4). In addition, from the high- 
energy perspective, the nn channel is (together with the 
three pion channel) essential to follow the spin in the 
Higgs di-tau decay channels at LHC 0, and thus to 
determine its spin and CP properties with the help of 
TAUOLA 0. 

2. Theoretical setting 

In t~ — > 7T~n°v T decays the electroweak part of 
the process is theoretically under control, while the 
hadronization of the quark currents is more involved 
since in the spanned energy region QCD is essentially 
non perturbative. The approximate chiral symmetry of 



light-flavoured QCD allows to build an effective field 
theory, known as Chiral Perturbation Theory OfPT) 
HI ED, which is able to describe accurately the 
low-energy part of the spectrum, but fails at larger in- 
variant masses 11 11 . In this region, new particles (p, 
K* , a,\, ...) are excited and their momenta and masses 
are large enough to prevent their use in the expansion 
parameters of the effective theory. Therefore these new 
active degrees of freedom have to be included in the ac- 
tion, and a new expansion parameter is needed. In this 
respect, the inverse of the number of colours of the QCD 
gauge group has proven to be a useful quantity to build 
the expansion upon fT2l[T3l[T4ll . A modelization of this 
idea in the meson sector for three flavours is provided by 
Resonance Chiral Theory (%T) H21 [161 [17] [HI , which 
recovers the^PT results at next- to-leading order (NLO). 
Since it is known that the lightest resonances dominate 
the dynamics, the infinite tower of states predicted in 
the Nc — > 00 limit can be restricted to the first excita- 
tions, taking into account as many states as required by 
the data. In addition, it is seen that the inclusion of res- 
onance widths is essential to describe the observed phe- 
nomenology, although widths arise at NLO in the 1 /Nc 
expansion. In principle, one can take into account this 
effect by computing off-shell widths consistently within 
R#T |fl9l [20l . On the other hand, as in any effective 
theory, the symmetry properties determine the operators 
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allowed in the Lagrangian, but leave the corresponding 
coupling constants unknown. However, the QCD short- 
distance behaviour of the Green functions and related 
form factors Q3] OH El EDO El |22] [23] HQ |25] Hi pro- 
vide a set of relations among these coefficients that ren- 
ders RyT more predictive. 

3. Vector form factor of the pion and fits to data 

Different approaches have been developed to deal 
with the diverse energy regimes. For s < M 2 , the com- 

[30] 



putation of F*(s) at NNLO in *PT El ISO 
proves useful. In order to enlarge the domain of appli- 
cability up to 1 GeV, unitarization techniques OTl [32l 
and the Omnes solution to the dispersion relation have 
been employed Il33l [34l [35l . Finally, in order to reach 
energies up to M T , the inclusion of the p' resonance 061 
and even a tower of resonances, inspired in the Nc — > °o 
limit 1 3711381 . have been proposed. From the RyT La- 
grangian, including only the p(770) multiplet one ob- 
tains 



Ml 



FUs) = 1 + 



F V G 



v^v 



F 2 



(2) 



where F is the pion decay constant in the chiral limit, 
My = M p , and Fy and Gy measure the strength of the 
pun and pV^ couplings, respectively, and V M stands for 
the quark vector current. If the vanishing of the form 
factor at large energies is required, the relation FyGy = 

F 2 is obtained, yielding Fy(s) = -j^rr- Now one can do 
better [33] and match this expression to the NLO result 
in^-PT. Final state interactions are included through the 
so-called chiral loop functions Ap(s,n 2 — M 2 ). Then, 
the unitarity and analyticity constraints determine the 
Omnes exponentiation of the full loop function, leading 
to 



M 2 y 
M 2 y-S 



exp 



96tt 2 F 2 



(3) 



Here one cannot simply include the resonance width by 
replacing M 2 ,- s by M v —s—iMyYy{s) in the propagator, 
since this would double count Bm[Ap(s)] and analytic- 
ity would be violated at NNLO in the chiral expansion. 
We follow instead a procedure similar to that proposed 
in Ref. |39 1 for the Kn vector form factor, in which uni- 
tarity and analyticity are satisfied exactly. The starting 
point is a form factor as in Eq. where the loop func- 
tions are resummed into the denominator: 



FUs) 



(0) 



My 



M v{ 1 + %^KW + 2 A ^)]}- S 



(4) 



M\ (l + ^rpi^e [A K ( S ) + |A*(i)]} - s - iMyYy(s) 

Thus the relevant {I - 1, J = 1) phase shift is taken to 
be 



tan<5 0) 



l 3mFJ(s) (0) 



%eFUs) {0) 



(6) 



This is now used as input for a three-subtracted disper- 
sion relation for the form factor. In this way one gets 

a 2 2 



F v (s) = expjais+ys 



s 3 r« d\(s') ) 

+ — ds' , 1 \. (7) 

n J4ml (s'y(s' - s-ie)) 

Kinematical isospin corrections can be easily in- 
cluded in Eq. (j7]l by considering different masses for 
the charged and neutral pions and kaons. In addition, at 
the same order one should also take into account elec- 
tromagnetic corrections |]2] |4D1 . which enter through a 
local term f£~°. and a global factor Gem(s) [41 1. Thus 
we consider three possible expressions for the form fac- 
tor to perform our fits to Belle data: 

• Fit I, corresponding to F y(s) from Eq. Q. 

• Fit II, same as I but with the inclusion of kinemat- 
ical corrections (m K ± + m n o, m.K* + m K o). 

• Fit III, including kinematical and electromagnetic 
corrections. 

Our fitting parameters are M p , F, a\ and It is found 
that without the inclusion of additional resonances, one 
can obtain good fits to Belle data for s < 1.5 GeV 2 . The 
best fit results to the first 30 points (central value of the 
bin corresponding to 1.525 GeV 2 , with 0.05 GeV 2 bin 
width) are shown in Table [T] where we have considered 
the l/N njI dN njI /ds distribution measured by Belle (this 
includes error correlations). These fits show, firstly, 
that the dispersive description of the form factor can 
indeed successfully account for the experimental data 
up to s < 1.5 GeV 2 , and secondly, that the approach 
employed by the Belle Collaboration (named here as 
II) is indeed an adequate one, as it yields the lowest 
X 2 Indf values according to our fits. Notice that, given 
the low energy threshold for this decay, the subtraction 
constants are fixed at a relatively low energy scale, and 
the dispersive representation turns out to be insensitive 
to the dynamics at large energies. In order to get a re- 
sult for the form factor that can be valid up to s = M 2 , 
the expression in Eq. |7]i can be e.g. matched at in- 
termediate energies to a phenomenologically adequate 
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function [ 1 1 like that given in Ref. [42] (included in the 
new version of TAUOLA [43 1), or the Gounaris-Sakurai 
parametrization |44] used by Belle [45|. 

In the table, the errors quoted in single brackets are 
those resulting only from the fit, i.e. neglecting the sys- 
tematic errors arising from our theoretical approach. 
These are e.g. given by the energy range to be fitted, the 
number of subtractions and the value of the upper inte- 
gration limit in the dispersive integral. In order to es- 
timate the associated total errors we have extended our 
fit up to energies in the range [1.325, 1.525] GeV 2 , we 
have taken into account the results for 2 and 4 subtrac- 
tions, and we have taken s M in the range [4, oo] GeV 2 . 
In this way we end up with the numbers quoted in dou- 
ble brackets in Table [T] Notice that the input values for 
the p mass and width still need to be translated to the 
physical pole values, which are reasonably lower [Q]. 

It is seen that our results show a lower x 1 l n d,f for 
our fitting options I and II, compared to that obtained 
for option III. Thus, the best agreement with the data is 
reached by including SU(2) isospin breaking only kine- 
matically, although comparable results are obtained in 
the isospin symmetric case. 

Finally, we have also computed the low-energy ob- 
servables (r 2 ^ and <f y appearing in the low-s expansion 
of Fy(s). Our results are quoted in Table |2j together 
with those obtained in previous analyses. It is seen that 
the values are entirely compatible, while the errors are 
found to be slightly reduced thanks to the good quality 
of present Belle data. 

4. Conclusions 

We have elaborated a dispersive representation of 
F y(s), which preserves analyticity and unitarity exactly 
and reproduces the low-energy limit of^PT up to NLO. 
We have performed different fits to Belle experimental 
data, which allow to determine our four input model 
parameters. The fits show a good agreement with the 
data, and no significant improvement is found after the 
inclusion of isospin breaking corrections. In addition, 
from our fits we have evaluated the low-energy quanti- 
ties (r 2 j and Cy, which turn out to be consistent with 
previous determinations. Our framework is also able to 
provide good quality fits to <x(e + e~ — > ti + jt~) scattering 
at low energies, which can be used to determine a™' L0 
from t decays and e + e~ scattering consistently. 
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Parameter 


Fit Value© 


Fit Value(II) 


Fit Value(III) 


M p 


0.8431(5)((17)) 


0.8280(4)((14)) 


0.8276(4)((21)) 


F n 


0.0901 (2)((5)) 


0.0902(2)((4)) 


0.0906(2)((4)) 


ct\ 


1.87(1)((3)) 


1.84(1)((3)) 


1.81(1)((2)) 


U2 


4.29(1)((7)) 


4.34(1)((7)) 


4.40(1)((6)) 


X'lndf 


1.37 


1.39 


1.56 




0.207(1)((3)) 


0.194(1)((3)) 


0.192(1)((4)) 



Table 1: Fit results to the Belle distribution, including correlations between errors. The errors in single and double 

brackets correspond to those arising only from the fit and those obtained after considering theoretical systematics, respectively. 
Energy units are given in GeV powers. T p (M p ) is obtained using the fitted values of M p and F n and is given only for reference. 



Determination 




4 (GeV- 4 ) 


Our fit 
Bijnens and Talavera [30] 
Pichand Portoles |34| 


10.86(14) 
11.22(41) 
11.04(30) 


3.84(3) 
3.85(60) 
3.79(4) 



Table 2: Low-energy observables of the vector pion form factor up to the quadratic term. We give the results from our fit, the 0(p 6 ) 
xVT analysis in Ref. 1 30 1 and the dispersive analysis in Ref. |34|. 



